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ABSTRACT 

The abundance of peaks in weak gravitational lensing maps is a potentially powerful cosmological tool, complementary to measure¬ 
ments of the shear power spectrum. We study peaks detected directly in shear maps, rather than convergence maps, an approach that 
has the advantage of working directly with the observable quantity, the galaxy ellipticity catalog. Using large numbers of numerical 
simulations to accurately predict the abundance of peaks and their covariance, we quantify the cosmological constraints attainable by 
a large-area survey similar to that expected from the Euclid mission, focusing on the density parameter, Dm, and on the power spec¬ 
trum normalization, erg, for illustration. We present a tomographic peak counting method that improves the conditional (marginal) 
constraints by a factor of 1.2 (2) over those from a two-dimensional (i.e., non-tomographic) peak-count analysis. We find that peak 
statistics provide constraints an order of magnitude less accurate than those from the cluster sample in the ideal situation of a perfectly 
known observable-mass relation; however, when the scaling relation is not known a priori, the shear-peak constraints are twice as 
strong and orthogonal to the cluster constraints, highlighting the value of using both clusters and shear-peak statistics. 
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1. Introduction 

Weak gravitational lensing (WL) is a powerful probe of large- 
scale structure, dark matter, and dark energy (e.g., Bartelmann 
&. Schneider 2001). In particular, cosmic shear surveys have 
demonstrated their ability to constrain cosmological parameters 
(e.g., Massey et al. 2007; Heymans et al. 2012), and the poten¬ 
tial of large-area shear surveys covering thousands of square de¬ 
grees to improve cosmological constraints to percent level ac¬ 
curacies is the prime motivation for ambitious programs like 
the on-going Dark Energy Survey (DES)*, and Stage IV ex¬ 
periments (Albrecht et al. 2006, 2009) like the Large Synoptic 
Survey Telescope^ (LSST, LSST Science Collaboration et al. 
2009), the Wide Eield Infrared Survey Telescope (WEIRST)^ 
(Spergel et al. 2013), and the Euclid^ (Laureijs et al. 2011) mis¬ 
sions. 

The shear correlation function (equivalently, power spec¬ 
trum), a measure of the second moment of the mass distribu¬ 
tion, is the standard tool for analyzing WL surveys. Other sta¬ 
tistical measures of shear maps incorporating higher order mo¬ 
ments are possible, and they become increasingly attractive in 
light of the significant gain in signal-to-noise expected from the 
planned large WL surveys. 


Send offprint requests to\ Nicolas Martinet, e-mail: martinet@iap. fr 
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In this paper, we consider the statistics of peaks in a shear 
map (or shear catalog) as a cosmological probe. We define shear 
peaks by filtering the map with an aperture designed to detect lo¬ 
calized projected mass concentrations, such as galaxy clusters. 
Indeed, WL surveys can be used to detect galaxy clusters, and 
the cluster counts then used as a cosmological probe (Kruse & 
Schneider 1999; Marian & Bernstein 2006). Projection effects, 
however, severely limit the purity of cluster samples defined 
through WL, despite attempts at constructing optimal filters, be¬ 
cause many peaks result from the alignment of small systems 
along the line of sight (White et al. 2002; Hamana et al. 2004; 
Hennawi & Spergel 2005). 

An alternative is to abandon the correspondence between 
shear peaks and clusters and simply use the statistics of peaks 
to characterize the projected mass distribution (Reblinsky et al. 
1999). This is the approach we adopt in the present work. These 
general shear peaks do not necessarily have any meaning as 
physical objects, being a combination of real clusters and chance 
alignments. Their abundance, however, like clusters, is sensitive 
to the underlying cosmology (Jain &. Van Waerbeke 2000; Wang 
et al. 2009; Dietrich & Hartlap 2010; Kratochvil et al. 2010). 

One disadvantage of this approach is that we do not posses 
a simple analytic form for the abundance of peaks as a func¬ 
tion of cosmological parameters. This is in contrast to the situ¬ 
ation with clusters, where practical expressions do exist for the 
mass function that greatly facilitate the theoretical prediction of 
cluster abundance and exploration of parameter space (Press & 
Schechter 1974; Jenkins et al. 2001 ; Tinker et al. 2008). We must 
instead resort to N-body simulations to predict WL peak abun- 
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dance, and we require large suites of simulations to explore the 
parameter space. 

Such studies have been undertaken by several authors in re¬ 
cent years (Dietrich & Hartlap 2010; &atochvil et al. 2010; 
Yang et al. 2011; Hilbert et al. 2012; Marian et al. 2012, 2013). 
In this paper, we perform a Fisher analysis of the constraints 
from peak counts in the context of upcoming Stage IV dark en¬ 
ergy surveys, comparing our results to constraints expected from 
cluster counts. We employ large suites of independent N-body 
simulations to mitigate what has been an important limitation of 
previous studies, and we work directly with shear measurements, 
rather than reconstructed convergence maps. 

The SUNGLASS code (Kiessling et al. 2011) is a rapid sim¬ 
ulation tool based on line-of-sight integration through N-body 
boxes to calculate the WL field. Its speed allows us to gener¬ 
ate large numbers of simulated shear maps, and hence determine 
peak abundance as a function of cosmological parameters and 
its variance. In particular, we are able to accurately calculate 
the derivative of peak abundance with respect to the parame¬ 
ters needed for the Fisher matrix. It is important to note that 
while many previous works have employed N-body simulations 
for similar analyses, they have relied on statistically shifted maps 
generated from a limited number of simulations (e.g., Dietrich & 
Hartlap 2010; Kratochvil et al. 2010; Yang et al. 2011; Hilbert 
et al. 2012; Marian et al. 2012, 2013). By contrast, the maps in 
this work are truly independent, with each map generated from 
a separate N-body realization. 

We work directly with ellipticity measurements, an unbiased 
estimator of the shear in the WL regime (e.g., Dietrich & Hartlap 
2010; Maturi et al. 2011; Hamana et al. 2012; Hilbert et al. 
2012), rather than convergence maps that have been used in sev¬ 
eral previous studies (e.g., Kratochvil et al. 2010; Yang et al. 
2011; Marian et al. 2012, 2013). Convergence is not the direct 
observable, but must be reconstructed from shear measurements. 
Our approach avoids the complexity added by this inversion. In 
addition, the use of shear greatly simplifies the nature of map 
noise, which is non-trivial to estimate in the case of the recon¬ 
structed convergence. 

Finally, we apply tomography to the peak statistics by di¬ 
viding the lensed background galaxies into redshift bins. This 
offers a two-fold advantage; first it allows us to remove fore¬ 
ground galaxies (z < 0.5) that tend to dilute the shear signal, al¬ 
though the mass distribution below this redshift still contributes 
to the statistics measured using only the higher redshift galaxies. 
Second, we can detect shear peaks in different redshift planes 
and examine the statistics both within and between planes. As 
with the shear correlation function, the additional radial infor¬ 
mation significantly increases precision on cosmological con¬ 
straints. Building on the work of Hennawi & Spergel (2005) and 
Dietrich & Hartlap (2010), we quantify the constraining power 
of tomographic shear peak statistics for Stage IV dark energy 
missions, such as Euclid. 

We begin by describing our WL simulations in Section 2. 
The peak detection procedure and its application to the sim¬ 
ulations are detailed in Section 3. In Section 4 we examine 
peak statistics and their use as a cosmological probe. Section 5 
extends the approach to tomography. We conclude with a fi¬ 
nal discussion and comment on future directions in Section 6. 
Throughout the paper, for concreteness, we consider the specific 
case of a survey similar to that of the Euclid mission with a fidu¬ 
cial flat ACDM cosmology specified by {Qm, h, crs, 

«s)={0.272, 0.728, 0.0449, 0.71, 0.809, 1.000) ( e.g., Hinshaw 
et al. 2013). 


2. Weak lensing simulations 

We employ the Fisher formalism to asses the cosmological con¬ 
straints expected from peak counts in Stage IV dark energy ex¬ 
periments, taking as typical characteristics those of the Euclid 
(Laureijs et al. 2011) mission. Euclid will survey 15,000 deg.^ in 
three infrared bands (Y, J, and H) and a single, wide optical fil¬ 
ter (combined riz bands), the latter with a point spread function 
(PSF) of 0.1 arcsec. The survey reaches a mean source galaxy 
density usable for lensing of 30 gals, arcmin^^ with a median 
redshift of z = 0.8, following the redshift distribution 

P(z)-(1) 

with zo - 0.7. 

Calculation of the Fisher matrix requires the derivative of 
the mean peak counts with respect to cosmological parameters, 
evaluated at the fiducial model. We also need the covariance of 
the peak counts about their mean in this model. Since we do 
not posses an analytical expression for the WL peak abundance, 
we must use simulations to calculate the expected peak counts. 
We need enough simulations to accurately determine the mean 
peak counts for each parameter variation and to determine their 
covariance in the fiducial model. This is non-trivial as much of 
the signal comes from non-linear scales in the mass distribution, 
which can only be properly modeled by N-body simulations. 

Simulation speed is therefore essential. In our study we use 
the SUNGLASS pipeline developed by Kiessling et al. (2011). 
We give a brief description of the SUNGLASS pipeline here but 
for a full prescription of how the SUNGLASS WL shear and con¬ 
vergence catalogs are generated, see Kiessling et al. (201 1). For 
a given set of parameters, SUNGLASS first generates an N-body 
realization with the GADGET-2 code Springel (2005). These sim¬ 
ulations are performed with 512^ paticles in a 512h '.Mpc box 
and the light cone is 100 square degrees and goes out to a redshift 
of z = 2.0. The convergence and shear of each mass point are cal¬ 
culated along the line of sight at multiple lensing source planes. 
This relies on the applicability of the Born approximation in the 
WL regime, i.e., that the mean path of the light bundle from 
a distant object remains adequately straight in the presence of 
lensing. The integration is much quicker than ray tracing through 
the simulation box, and allows the production of large suites of 
simulated WL observations. The shear and convergence is then 
interpolated on to the individual particles in the light cone, pro¬ 
viding a highly sampled catalog of shear and convergence along 
the line of sight. 

The speed of SUNGLASS enables us to produce many in¬ 
dependent realizations, an improvement over previous studies 
that had to resort to shuffling the results from single realiza¬ 
tions. The final WL catalogs are constructed by down-sampling 
the highly sampled SUNGLASS shear and convergence catalogs 
to 30 galaxies per square arc minute using the source galaxy 
redshift distribution of Eq. (1), and assuming that the galaxies 
trace the dark matter exactly. To model the shape noise arising 
from intrinsic galaxy ellipticity, we add a random ellipticity to 
each source galaxy in the simulated shear catalogs according to 
a Gaussian distribution with zero mean and standard deviation 
cr^. In practice, we use a Box-Muller Gaussian random number 
generator, drawing two uniform numbers (x, y) over the interval 
[0,1) to obtain two Gaussian random numbers (ei, € 2 ) that we 
add, respectively, to the shear components, yi and j 2 - 
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ei - sJ-2 In (x) cos (27ry) 

V2 


( 2 ) 
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62 = V“2 In (x) sin (2;ry). (3) 

V2 

Null values in the first step are replaced by arbitrarily chosen 
values. 

This procedure corresponds to drawing a tangential elliptic- 
ity from a Gaussian of zero mean and dispersion cr^, which we 
take to be cr^ = 0.3 (e.g., Leauthaud et al. 2007). We note that 
the shear components, yi and 72 , are typically only a small per¬ 
cent of the noise, ei and 62 - Under the assumption that the noise 
is uncorrelated (i.e., in the absence of important intrinsic align¬ 
ment effects), the mean ellipticity averaged over a number of 
source galaxies approaches the shear signal. We will discuss this 
fundamental hypothesis in Section 5. 


3. Peak detection 

3.1. Method 

To identify peaks in a shear catalog, we employ aperture mass 
filtering (Schneider et al. 1998; Bartelmann & Schneider 2001) 
with an outer annulus to remove the integration constant, thereby 
resolving the finite space inversion problem. This technique does 
not return the true mass within the aperture, for example, when 
centered on a cluster; for that, one would have to know the 
true shape of the mass distribution as seen through the filter. By 
adopting a mass distribution, such as the NFW profile (Navarro 
et al. 1997), one can determine the true mass and also find clus¬ 
ters in WL surveys (e.g., Marian et al. 2012, 2013). Our goal, 
however, is not to measure the true mass of physical objects, but 
to compare the number of peaks expected in a WL survey for 
different cosmologies. 

The aperture mass can be calculated either from the conver¬ 
gence field, K, or the shear field, 7 . In the former case, the aper¬ 
ture mass is calculated by integrating the convergence within the 
aperture centered at position Oq (a two-dimensional vector in the 
plane of the sky), 

Map(0o) = J d^o U{\e - eoMe), (4) 

where t/ is a filter chosen to best fit the lens mass density profile. 

In a WL survey, the observable is actually the source galaxy 
ellipticity, e. It is related to the reduced shear, g, by 

<f> = ^ = 7 -^ ~ r> (5) 

which tends to the shear, 7 , in the weak lensing regime where 
{k,j) « (1,1). The indicated average is over random orien¬ 
tations of intrinsic galaxy ellipticity. A WL survey thus directly 
measures shear, 7 , rather than the convergence. Working directly 
with the observable quantity, 7 , we avoid the non-trivial step of 
integrating the shear to obtain the convergence field, which is a 
derived quantity. 

Shear and convergence are two mathematically distinct, al¬ 
though related, quantities, the first being a spinor of spin two 
and the second a scalar on the sphere. To adapt the expression 
for the aperture mass to the case of shear, we first define the 
scalar tangential shear, 7 ,, for a galaxy image that accounts for 
both components of the shear (71 and 72 ), 

7 t( 0 , 60 ) = - [yi{ 0 ) cos 2 (^( 0 , 0 o) + JiiO) sin 2 (j){ 0 , 0 o)], ( 6 ) 

where is the angle giving the position of the galaxy image { 0 ) 
relative to an arbitrary fixed axis running through the center of 
the aperture, at position 0 q, in the image plane; this fixed axis 


defines a local cartesian coordinate system in the plane of the 
sky with origin positioned on the filter center. The aperture mass 
equation can then be rewritten in terms of the tangential shear 
and the new filter function Q: 

2 r® 

^ cfG'e' u{\ff\)-u{\0\), (7) 

Map(0o) = Juf-ee(|e-0ol)rt(0,0o). (8) 

The convergence has been used in several previous weak 
lensing peak studies (e.g., Kratochvil et al. 2010; Yang et al. 
2011; Marian et al. 2012, 2013), although recent studies increas¬ 
ingly work directly with the shear (e.g., Dietrich & Hartlap 2010; 
Maturi et al. 2011; Hamana et al. 2012; Hilbert et al. 2012), 
which is also the approach adopted in this paper. 

We use the following weight functions, U for conver¬ 
gence and Q for shear, appropriate for a circular aperture (see 
Bartelmann & Schneider 2001): 


t/(0) = 


0 


ifO < 6 > < 0 y 
6 \ < 6 < 02 
elsewhere 


Q(0)^ 


0 


if til < 0 < 02 
elsewhere 


(9) 

( 10 ) 


In practice, the integral over the aperture weight function be¬ 
comes a sum weighted by the number density of galaxy images, 
n. 


Map(eo) = - y Ui\0i - 0o\M0i) (11) 

n ^ 

I 

Map(0o) = - y Q(\0i - eo|)7t(0/, 0o), (12) 

n ^ 

I 

where 7 /( 0 ,, 0 o) and K{0i) are the tangential shear and the conver¬ 
gence of the image at 0 , relative to the point 0 o. 

By expressing the equations in terms of discrete sums, we 
explicitly account for the sampling inherent in the observations, 
i.e., we only measure the shear where there is a source galaxy. 
Moreover, we use the actual number density of galaxies in the 
aperture, rather than a fixed, average value. The link from simu¬ 
lations to observations is trivial, as it is sufficient to replace the 
tangential shear, 7 ,, in the last equation by the tangential com¬ 
ponent of the ellipticity, e,. This is the principal interest of using 
shear instead of convergence. 

The aperture mass is convenient because in the case of the 
shear, it permits simple calculation of its variance due to the in¬ 
trinsic ellipticity of the source galaxies, cr^: 


o-(A7ap) = 


V2n 


xl/2 


V i y 


(13) 


This allows us to define a local noise level and peak detection 
threshold, which is another strong argument in favor of using 
shear peaks rather than convergence. We then define peak am¬ 
plitude as 


^ap ^ ZiQ(\Oi-Oo\)7,iOi,Oo) 


(14) 
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The amplitude of a peak is independent of the normaliza¬ 
tion of the weight function, but does depend on the number of 
galaxies in the aperture through both the sums in numerator and 
denominator. Peak amplitude thus varies as the square root of 
the number of galaxies in the aperture. Intrinsic galaxy ellip- 
ticity also affects the signal-to-noise value, although it does not 
change the relative intensity of the peaks. As mentioned, we take 
cTf = 0.3 (Leauthaud et al. 2007). 

A critical point is the size and shape of the aperture. The 
chosen shape will preferentially select a specific form of struc¬ 
ture, such as clusters or filaments, while the size will favor one 
cluster size over others. We must also adapt the aperture to in¬ 
clude enough galaxies to optimize the signal-to-noise over the 
random shape noise (this point will be discussed in the section 
on tomography). 

For simplicity, we adopt a radially symmetric aperture of 
fixed angular radius. The inner radius is set to 9i - 3.5', corre¬ 
sponding to the typical size of a cluster at redshift z - 0.3, where 
the contribution to shear peaks is the most important (Dietrich & 
Hartlap 2010), to exclude any contamination by cluster galaxies 
and the strong lensing regime, while the outer radius is set to 
02 = 10', which is roughly the limit of the lensing effect at this 
redshift (see Hamana et al. 2012). In practice, we have found that 
the mass inside a given aperture strongly depends on the inner 
radius. In a future work, we plan to use a set of aperture sizes 
to extract information from different scales. One could also use 
an adaptive matched filter to preferentially select galaxy clus¬ 
ters (e.g., Marian et al. 2012). This is not our goal in this first 
study, and we leave the identification of an optimal filter to a fu¬ 
ture work. Finally, peaks are selected to be larger than all their 
neighbors in a radius equal to that of the aperture. Peaks situated 
at less than 62 from the map edges are discarded as they are not 
computed in the proper aperture. 

4. Peak statistics 

We first present our statistical methodology and results from a 
non-tomographic analysis of the peak counts. Section 5 then ex¬ 
tends the analysis to tomographic peak counts. 

4.1. Method 

We implemented two statistical measures: a test and the 
Fisher information matrix, both defined over bins of peak height. 
We chose our bins to include the same number of peaks based on 
the mean peak counts in our fiducial cosmology. This bin size is 
then maintained for the other cosmologies. Bin widths are given 
in Table 1, along with the number of peaks for one fiducial real¬ 
ization. 

Let Ni,r be the number of peaks in bin i of realization r for 
a given cosmology, and R be the total number of realizations of 
this cosmology. Defining (A), as the mean number of peaks in 
bin i, averaged over all R realizations, we calculate the covari¬ 
ance matrix of the binned peak counts as 

r - <N}i)(Nj,r - (A),) if / ^ 

I Zf (Nu - {N}i){Nj,, - {N}j) if / = ;• ’ ^ ^ 

using R - 150 independent realizations of the fiducial cosmol¬ 
ogy. Each of these realizations corresponds to a lightcone of 
lOOdeg^, and we subsequently normalize the covariance ma¬ 
trix to an area of 15,000deg“ (e.g., the useable extragalactic 
sky and Euclid target). Figure 1 shows the correlation matrix. 


Table 1. Bin widths and peak counts for one fiducial realization. 
The second column gives the lower bound on the peak signal-to- 
noise for each bin and the third the total number of peaks in the 
bin; for example, the first line reads: in bin 1 there are 70 peaks 
with signal-to-noise between 3.0 and 3.2. 


bin number 

bin size 

N 

1 

3.0 

70 

2 

3.2 

72 

3 

3.4 

65 

4 

3.5 

67 

5 

3.7 

72 

6 

3.9 

74 

7 

4.2 

69 

8 

4.5 

77 

9 

5.0 

66 

10 

5.7 

90 


i.e., the covariance matrix normalized to unity along the diago¬ 
nal. We see that the correlation between bins is less than 20% 
except for the higher signal-to-noise bins where it can reach 
up to 40%. It seems reasonable that the stronger peaks would 
be more correlated between bins, with signal being dominated 
and seen by successive source planes, while the lower signal- 
to-noise peaks would be more affected by noise variations and 
projections between source planes. This agrees with the fact that 
high signal-to-noise peaks mostly correspond to galaxy clusters, 
while low signal-to-noise peaks are dominated by projections of 
large-scale structure and noise, as shown by Maturi et al. (201 1). 



Fig. 1. Correlation matrix of the shear peak distribution for a 100 
square-degree field, calculated by averaging over 150 realiza¬ 
tions of the fiducial cosmology (see Eq. 15). This is the covari¬ 
ance matrix normalized to unit diagonal. The x and y-axis values 
are the peak height bin numbers, defined in Table 1 . 

This covariance matrix is then used to compute either the 
(Eq. 16) or the Eisher matrix (Eq. 17). Thex^ is expressed as 

Ar, f) = ^ (A{; - {Nf)d{Clr\N{^ - {Nf)j), (16) 

U 

where / represents the fiducial cosmology and /' a reference 
cosmology. We note that /' can be the same as / if we wish 
to compare one realization of a cosmology to all the other real¬ 
izations of the same cosmology. This allows us, in particular, to 


4 










Martinet et al.: Constraining cosmology with shear peak statistics 


test whether the^^ variable is actually distributed according to a 
;^f^-distribution. The Fisher matrix is given by 


ft.Pb 


^ d/7a apb 


(17) 


where and p\, are two cosmological parameters. We note 
that this expression is easily interpretable only in the case of 
Gaussian distributed bin counts, i.e., as the number of peaks be¬ 
comes large. Cosmological parameter constraints are then ob¬ 
tained by inverting the Fisher matrix, 

(18) 


Following Hartlap et al. (2007), we use the unbiased estima¬ 
tor, C ', for the inverse covariance matrix. Under the assumption 
of Gaussian errors and independent data vectors, it is related to 
the inverse of the estimated covariance matrix (Eq. 15) through 
the number of realizations, R, and the number of degrees-of- 
freedom, D, as 


c-i = 


R-D-2 

R-1 




(19) 


where in our case D is the number of peak-height bins. This cor¬ 
rection will thus be more important in our tomographic analy¬ 
sis where we build the covariance matrix through the assembled 
peak distributions of several redshift slices. 

The derivative of the peak counts with respect to cosmologi¬ 
cal parameters averaged over all the realizations r is given by 


d{N)i _ 1 y NiAp^ + Apa) - NiApa - Apa) 
9Pa R ^ 2Apa 


( 20 ) 


where Apa is the variation of the cosmological parameter pa- 
This calculation requires a sufficient number, R, of numeri¬ 
cal simulations of each cosmology to accurately determine the 
derivatives. For a single parameter variation, we used 250 cos¬ 
mological simulations (150 for the fiducial cosmology and 50 
for a variation of -i-Apa and 50 for a variation of -Apa). 

To determine if this is sufficient, we perform a convergence 
test on numbers of realizations by comparing the derivatives of 
the peak counts for increasing numbers of realizations of the 
modified cosmologies. In Fig. 2 we show the derivatives with 
respect to and erg when varying the number of realizations 
from 10 to 50 in increments of 10 realizations. We note that the 
derivatives do not significantly evolve beyond 30 realizations, 
justifying our choice of /? = 50 for the modified cosmologies. 

We use a larger number of realizations of the fiducial cosmol¬ 
ogy because the covariance matrix is computed for that model, 
while only the mean peak counts are required for the other cos¬ 
mologies. Taylor & Joachimi (2014) calculated the accuracy of 
the covariance matrix given the number of realizations and of 
degrees of freedom of the data vectors. Following their Eq. (13), 
we estimate the precision of our covariance matrix to be better 
than ~ 13% with our 150 realizations. The additive factor of 
2v“^, where v is the desired accuracy, in the required number 
of realizations limits in practice the achievable accuracy on the 
covariance matrix; sub-percent accuracy, for instance, would de¬ 
mand at least 40,000 realizations. Our choice of 150 seems rea¬ 
sonable for the present test-study, but this issue calls for further 
attention and presents a crucial difficulty for many dark energy 
probes based on large-scale structure, such as cosmic shear. 



Fig. 2. Derivatives of the peak counts with respect to Om (top) 
and cTg (bottom) as defined in Eq. (20) as a function of the 
number of realizations of the modified cosmologies. Results 
for R = (10,20,30,40,50) correspond, respectively, to the red, 
green, blue, purple and cyan curves. There is little change be¬ 
yond R = 30. 


4.2. Results 

As an illustration of peak count statistics, we study achievable 
constraints on two cosmological parameters: the total matter 
density, and the present-day linear matter power spectrum 
normalization, erg. In the standard ACDM model, these are well 
constrained by cosmic microwave background (CMB) observa¬ 
tions (Hinshaw et al. 2013; Planck Collaboration et al. 2013a). 
Methods measuring their values at low redshifts, such as peak 
counts or other gravitational lensing observations, are useful to 
search for extensions of this simple model. The power spectrum 
normalization, erg, is a good example: differences between val¬ 
ues obtained from the CMB and low redshift methods could in¬ 
dicate the need for a non-minimal neutrino mass (e.g., Planck 
Collaboration et al. 2013b; Rozo et al. 2013; Battye & Moss 
2014). 

Each parameter was varied by 10% from its fiducial value 
given in Table 2 to define a reference cosmology. We generated 
150 realizations of the fiducial cosmology and 50 realizations 
of each reference cosmology. When varying one parameter, the 
other remains at its fiducial value. However, when varying the 
matter density parameter, Q^, the dark energy density parame¬ 
ter, Da, was also adjusted in order to maintain a fiat Universe, 
Gat + Ga = 1. 
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Table 2. Cosmological parameter values for the fiducial and ref¬ 
erence cosmologies. The quantity R is the number of realizations 
of each cosmology. When varying one parameter, the other re¬ 
mains fixed at its fiducial value. 



Hm 

cm 

R 

fiducial 

0.272 

0.809 

150 

low 

0.245 

0.728 

50 

high 

0.299 

0.890 

50 


4.2.1. chi-squared distribution 

We first examine the distribution of values (Eq. 16) in the 
fiducial cosmology, using the covariance matrix calculated over 
the 150 realizations and dividing the peak heights into 10 bins 
of equal numbers of peaks, on average. The top panel of Fig. 3 
shows that the distribution observed in the simulations is reason¬ 
ably well represented by a true;^^^ distribution with 10 degrees of 
freedom, although with a slight deviation manifest by the some¬ 
what larger variance. 

When comparing a modified cosmology to the covariance 
matrix of the fiducial cosmology, we see that the distribution 
strongly diverges from a true;t'“ law (lower panel of Fig. 3). This 
result illustrates the potential of this method to constrain cosmo¬ 
logical parameters. The next step is to compute the constraints 
we would achieve with the Euclid survey using the Fisher for¬ 
malism. 

4.2.2. Fisher information 

Following Eq. (17) we compute the Fisher matrix and invert 
it to obtain constraints on the cosmological parameters. Two- 
dimensional constraints are plotted in Fig. 4 with Icr and 2cr con¬ 
fidence contours. These constraints are summarized in Table 4. 


5. Tomography 

We develop a tomographic approach to peak count statistics 
with the aim of exploiting the radial information by dividing the 
source galaxies into redshift bins and detecting peaks to each 
source plane separately. We then perform a joint statistical anal¬ 
ysis of the multi-plane peak counts. For example, with two red- 
shift bins we would have a peak distribution consisting of 20 
bins in which the first 10 bins represent the distribution of peaks 
detected to the first source plane, and the remaining 10 the dis¬ 
tribution to the second source plane. Our analysis employs the 
full covariance of these 20 bins. We note that this approach dif¬ 
fers from that employed by Dietrich & Hartlap (2010) in that we 
do not attempt to localize individual peaks in redshift space; the 
two approaches, however, access the same information. An im¬ 
portant issue with tomography is to ensure that we have enough 
galaxies in each aperture for the average ellipticity noise to be 
negligible compared to the average tangential shear. 


5. 1. How to slice the redshift dimension 


Adopting a Gaussian random distribution of intrinsic galaxy el- 
lipticities with zero mean and dispersion cr^, the shape noise over 
an aperture is 


rr ap 


o-g 

V2^’ 


( 21 ) 


where Aap is the average number of galaxies in the aperture. 
We denote y as the desired ratio between the average tangen- 



0 5 10 15 20 25 30 35 40 45 

Fig. 3. Normalized x^ distribution based on the covariance ma¬ 
trix from the 150 realizations of the fiducial cosmology. We con¬ 
sider 10 bins in peak height with equal numbers of peaks in each 
bin (on average). Black histograms represent our data and red 
curve represents a theoretical distribution with 10 degrees of 
freedom. The top panel shows the distribution for 150 fiducial 
realizations. The bottom panel is the distribution for 50 realiza¬ 
tions with Qn, increased by 10% from its fiducial value. 



Fig. 4. Fisher joint conditional constraints on Qm and erg with 
non-tomographic peak count statistics for a Euclid-like survey; 
the red and black contours delineate 1- and 2-cr significance lim¬ 
its. 


tial shear, ('yt)ap, and the aperture shape noise. An estimate of 
the number of source galaxies required per aperture is then 


Nap 


yo-e 


K V2('yt)ap J 


( 22 ) 
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We take a shape-noise dispersion of CTf = 0.3 and an average 
shear value of 0.04 (e.g., Jain & Taylor 2003; Jain & Seljak 
1997). These values with y -1 yield a required number of about 
1400 source galaxies per aperture. 

Working with fixed aperture size, the available number of 
galaxies per aperture depends on redshift. Using the distribution 
of redshifts in our simulations (Eq. 1) normalized to the mean 
Euclid galaxy density of 30 galaxies per square arc-minute, we 
can estimate the number of galaxies per aperture for any slice of 
redshift. For the most distant redshift bin, the condition of having 
at least 1400 galaxies per aperture is satisfied for 1.43 < z <2. 
We adopt the same number of galaxies per redshift slice to avoid 
favoring any particular redshift bin. The condition on a mini¬ 
mal number of galaxies per aperture then directly translates into 
a condition on the maximum number of redshift bins. We also 
do not use the z < 0.5 redshift range to avoid diluting the sig¬ 
nal. These conditions allow us to perform a tomographic analysis 
with up to five redshift slices between redshift 0.5 and 2. We note 
that relaxing the constraint on the shear to ellipticity ratio would 
allow more redshift slices. We also note that this approach is lim¬ 
ited by the uncertainty on the photometric redshift information, 
which is on the order of cr(z) = 0.05 x (1 + z). 

For a first tomographic study, we use the following five red¬ 
shift slices with equal numbers of galaxies: 0.5 < z < 0.73, 
0.73 < z < 0.93, 0.93 < z < 1.15, 1.15 < z < 1.43, and 
1.43 < z < 2. The mean density in each slice is about five galax¬ 
ies per square arc-minute. This corresponds to about 1400 galax¬ 
ies per aperture and a shear to ellipticity ratio of about seven. 

5.2. Results 

We use the same simulations and shape noise realizations as in 
Section 4 when studying the two-dimensional peak counts. The 
size of the peak-amplitude bins is determined to include the same 
number of peaks in each bin of a given redshift slice, based on 
the mean peak counts in our fiducial cosmology. Fig. 5 shows the 
full correlation matrix across all bins and source redshift planes, 
with the first 10 bins corresponding to the lowest redshift source 
plane and followed in sequence out to the highest redshift plane. 

Table 3 gives the mean and one sigma variation of the num¬ 
ber of peaks over the 150 realizations of the fiducial cosmology. 
The mean number of peaks are also shown in Fig. 6. 

As seen from Table 3 and Fig. 6, peaks detected toward a low 
redshift source plane tend to also be detected when using higher 
redshift planes. This correlates the peak counts between source 
planes, as indicated by the non-zero off-diagonal elements of the 
correlation matrix, especially in the higher signal-to-noise bins. 
The bins are not fully correlated, however, because new peaks 
are detected beyond the lower redshift source planes as we move 
outward. This tomographic view of the peak distribution con¬ 
tains valuable cosmological information that increases the con¬ 
straining power of the peak counts. 

This can be appreciated from the differences in the distri¬ 
butions shown in Fig. 7. The black histogram in the upper panel 
gives the distribution in the simulations, compared to a true x^ 
distribution with 50 degrees-of-freedom traced by the solid red 
line. The lower panel gives the distribution of our;^^^ variable for 
the same non-fiducial cosmology considered in the lower panel 
of Fig. 3. As before, the observed histogram strongly differs from 
the pmex" distribution. The fact that the histograms differ even 
more than in the two-dimensional case illustrates our increased 
ability to distinguish these cosmological models. 

We quantify this gain with a Fisher analysis of the same cos¬ 
mological parameters considered in the two-dimensional case. 



5 10 15 20 25 30 35 40 45 50 


Fig. 5. Correlation matrix of the shear peak distribution for a 
100 square-degree field with tomography, calculated by averag¬ 
ing over 150 realizations of the fiducial cosmology (see Eq. 15). 
Each redshift slice is divided into ten bins of peak height. This is 
the covariance matrix normalized to unit diagonal. Correlations 
between bins are less than 20%, except at the highest signal-to- 
noise. 



S/N 


Fig. 6. Mean number of peaks over the 150 realizations of the 
fiducial cosmology. The red histogram corresponds to peaks de¬ 
tected in the 2D analysis while green, blue, pink, cyan, and yel¬ 
low, respectively, correspond to peaks detected in the 0.5 < z < 
0.73, 0.73 < z < 0.93, 0.93 < z < 1.15, 1.15 < z < 1.43, and 
1.43 < z < 2 redshift slices. 


The increase in the number of effective degrees-of-freedom in 
the peak distribution drops the precision on our estimated covari¬ 
ance matrix to 15%, calculated according to Taylor & Joachimi 
(2014). Doubling the number of realizations to R - 300 would 
provide 10% precision, a small gain compared to the computa¬ 
tional time required to generate twice as many realizations. We 
also verify that R - 50 realizations of the modified cosmologies 
is sufficient for calculation of the derivatives of the mean counts: 
as before, they are stable beyond 30 realizations. The constraints 
from our tomographic analysis are given in Fig. 8 and listed in 
Table 4, and they are compared to the two-dimensional case in 
Fig. 9. 
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Table 3. Mean and one sigma variation of the number of peaks over the 150 realizations of the fiducial cosmology. 



3<^<5 

5<^<7 

7<^<9 

9 < - 

3 < - 

All galaxies (0.5 < z < 2) 

572 ± 19 

104 ± 12 

18 ±5 

7±3 

702 ± 20 

0.5 < z < 0.73 

344 ± 16 

3±2 

0±0 

0±0 

347 ± 16 

0.73 < z < 0.93 

392 ± 16 

7±3 

0± 1 

0±0 

399 ± 16 

0.93 <z< 1.15 

445 ± 18 

12±4 

1± 1 

0±0 

458 ±19 

1.15 <z< 1.43 

494 ± 19 

17 ±5 

1 ± 1 

0±0 

512 ±20 

1.43 <z<2 

553± 16 

27 ± 6 

2± 1 

0±0 

583 ± 17 





Fig. 7. Normalized distributions with respect to the covari¬ 
ance matrix for 150 realizations of the fiducial cosmology with 
five source redshift planes. Each of the five source planes is as¬ 
sociated with 10 signal-to-noise bins containing equal numbers 
of peaks. Black histograms trace the distribution observed in the 
simulations and the red curve represents a theoretical x^ distri¬ 
bution with 50 degrees-of-freedom. The top panel gives the dis¬ 
tribution for the 150 realizations of the fiducial model. The lower 
panel gives the distribution for the 50 realizations of the alternate 
cosmology with Qm increased by 10% from its fiducial value. 

6. Discussion 

A number of authors have investigated the potential of lensing 
peak counts as a cosmological probe, although few have consid¬ 
ered large Stage IV projects such as Euclid (Yang et al. 2011; 
Maturi et al. 201 1; Hilbert et al. 2012; Marian et al. 2012, 2013). 
Kratochvil et al. (2010), for example, examined the difference in 
the x^ distributions from different cosmological models; how¬ 
ever, it is difficult to make a direct comparison with our re¬ 
sults since they varied a different set of cosmological parameters. 
In their pioneering study of tomographic peak counts, Dietrich 
& Hartlap (2010) considered a CFHTLS-like survey of about 
180deg2. 


Fig. 8. Predicted joint conditional constraints on Qm and erg for a 
Euclid-\\]/x survey using tomographic peak count statistics. Solid 
red and dashed black contours correspond to Icr and 2cr confi¬ 
dence regions, respectively. 

Table 4. Predicted cosmological parameter constraints for a 
EuclidAikR survey. The numbers give Icr uncertainties and the 
corresponding relative percentages of the fiducial parameter val¬ 
ues Qn, = 0.272 and erg = 0.809. 




<5o-8 

Unmarginalized 

All galaxies 

0.0012 (0.43%) 

0.0018 (0.22%) 

Tomography 

0.0010 (0.35%) 

0.0014(0.17%) 

Marginalized 

All galaxies 

0.0037 (1.34%) 

0.0056 (0.69%) 

Tomography 

0.0018 (0.66%) 

0.0026 (0.32%) 


In this work we focus on a typical Stage IV survey charac¬ 
terized by the redshift distribution of a Euclid-\]k£ survey, using 
a suite of independent numerical simulations to examine pos¬ 
sible cosmological constraints. Figure 9 and Table 4 quantify 
gains in constraining power by using tomography. We improve 
the marginal constraints on Qm and erg by more than a factor of 
two over the two-dimensional (non-tomographic) analysis. As to 
be expected, the conditional constraints are improved by smaller 
factors; about 1.2 for Qm and 1.3 for erg. 

Among previous studies of peak-statistics, only Dietrich & 
Hartlap (2010) and Yang et al. (201 1) have applied tomography, 
the former maximizing the peak signal-to-noise given the red¬ 
shift distribution of galaxies, and the latter placing source galax¬ 
ies at either = 1 or = 2 and ray-tracing through simulations. 
Dietrich & Hartlap (2010) demonstrated the gain from tomo¬ 
graphic peak counts for a survey of 180deg“. In their FSST-like 
survey of 20,000deg^, Yang et al. (201 1) noted an improvement 
of the marginal constraints by factors of two to three when using 
tomography, in qualitative agreement with our results. 
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0.265 0.27 ^ 0.275 ^ 0.28 


Fig. 9. Fisher ellipses (at Icr ) for and erg for a Euclid-Yik& 
survey. The blue dashed curve shows the joint conditional con¬ 
straints without tomography, while the solid red contour gives 
those when using tomography with five redshift bins, demon¬ 
strating the important gain. 


We also note that the impact of shape noise is effectively re¬ 
duced with tomography, in particular for the peaks generated by 
structures at the higher redshifts. Binning the sources into red- 
shift planes removes the shape noise contributed by foreground 
galaxies that do not carry any signal on the higher redshift peaks. 

We compare our conditional constraints to those obtained by 
other authors. These vary over the range 0.0006 < bn„ < 0.0009, 
according to Hilbert et al. (2012) and Marian et al. (2012, 2013), 
and 0.0013 < bo-g < 0.0016, according to the same authors 
and Maturi et al. (2011). While these studies differ in a num¬ 
ber of respects, the agreement on the conditional constraints 
among these authors and our results is very good. Indeed, we 
reach (bn,„,bo-g) = (0.0012,0.0018) without tomography, and 
(0.0010,0.0014) with tomography. The very small difference 
from the literature can be attributed to our use of a slightly lower 
survey area: 15,000deg^ compared to 18,000 to 20,000deg“ in 
these other studies. 

Finally, we note that some authors have examined constraints 
on other cosmological parameters. In particular, it has been 
found that shear peaks have a good ability to constrain the dark 
energy equation-of-state w (e.g., Yang et al. 2011; Hilbert et al. 
2012; Marian et al. 2012, 2013) and primordial non-Gaussianity 
/nl (e-g-. Maturi et al. 2011; Hilbert et al. 2012). 

We also compare our forecasted constraints with those from 
other cosmological probes, in particular from cluster abundance 
studies. While some lensing peaks do arise from individual 
clusters, peak statistics represent a more general description of 
the matter distribution because many originate from projections 
along the line of sight. Figure 10 compares our predicted con¬ 
straints from tomographic peak statistics with those from current 
galaxy cluster constraints as summarized by Allen et al. (2011). 
The peak counts yield much tighter constraints than those from 
the current cluster analyses, which is not surprising because we 
are comparing present day cluster constraints to future lensing 
counts. 

A more pertinent comparison is between constraints pre¬ 
dicted from cluster photometric sample for the Euclid survey 
and the peak count constraints. This is shown in Fig. 11. The 
cluster constraints have been evaluated by Sartoris et al. (2015) 
for a Euclid cluster catalog considering the information provided 
by cluster number counts. We note that, unlike what has been 
done in Sartoris et al. (2015), constraints from clusters have 


been performed by varying only the erg and cosmological 
parameters and the 4 parameters that describe the bias, the scat¬ 
ter of the observable mass relation, and their redshift evolution. 
This has been done to compare in a more appropriate way the 
constraints obtained from the clusters and those from the shear 
peaks. The reduction in the number of free cosmological param¬ 
eters explains why the cluster constraints shown in Fig. 11 are 
smaller than those shown in Sartoris et al. (2015). We see that 
the constraints from tomographic peak counts are weaker than 
cluster counts by an order of magnitude when supposing that 
the observable-mass relation is fully known a priori (blue dotted 
ellipse). However, the shear-peak constraints are almost twice 
as strong when not making any assumption on the scaling rela¬ 
tion parameters and their evolution (green dash dotted ellipse). 
In addition, it is worth noting that the shear-peak contours are or¬ 
thogonal to the clusters when the observable-mass relation is not 
known a priori. This essentially shows the value of using both 
clusters and shear-peak statistics. As it is difficult to predict how 
well we will be able to constrain the scaling relation, we show 
here two extreme cases for the cluster constraints, with the idea 
that the observational constraints should lie somewhere between 
the blue and green ellipses. In particular, the green contours are 
very conservative, as we could in principle already constrain the 
scaling relation at redshift z = 0, which would reduce the errors 
on the cosmological parameters. 

We expect peak statistics to complement more standard two- 
point lensing statistics. In fact, the study by Dietrich & Hartlap 
(2010) suggests that peak statistics could yield tighter constraints 
than the classical 2-point correlation function, a result expected 
given that the peak statistics contain higher order correlations, 
and later confirmed by, e.g., Marian et al. (2012, 2013). Peak 
counts would also appear to be less affected by shape measure¬ 
ment systematics than the shear power spectrum in the sense that 
it is more difficult to reproduce the pattern necessary for peak 
identification than to affect the amplitude of two-point correla¬ 
tions. In similar vein, we expect that their respective sensitivi¬ 
ties to photometric redshift uncertainties will not be the same. 
Overall, the two methods will not share the same systematics 
and therefore offer important complementarity. 

7. Conclusion 

We have found that shear peak statistics offer a potentially pow¬ 
erful cosmological probe, in agreement with previous studies. 
As an advance along these lines, our results clearly illustrate the 
gain of using tomography in the framework of Stage IV dark 
energy surveys, i.e., separating the source galaxies into redshift 
planes and counting peaks to each plane. With tomography, we 
improve the conditional (respectively, marginal) constraints by a 
factor of 1.2 (resp. 2) on Qm and erg. 

For a large-area survey, typified here as that from a Euclid- 
like mission, we estimate that the peak-count constraints are an 
order of magnitude less powerful than those predicated from 
galaxy cluster evolution when the observable-mass relation is 
fully known a priori, while they are twice as strong when not 
making any assumption on this relation. The peak counts, how¬ 
ever, have the great advantage of not relying on such a scaling 
relation that may prove difficult to establish to high accuracy. 

We have only explored the two parameters Qm and erg in 
the present study, but plan to extend to other parameters, in¬ 
cluding the dark energy equation-of-state and primordial non- 
Gaussianity in future work. Further topics warranting explo¬ 
ration include the impact of various systematics, such as intrinsic 
alignments, photometric redshift errors, and shape measurement 
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Fig. 10. Constraints from tomographic peak counts compared 
to current galaxy cluster constraints, with contours giving the 
1- and 2-cr confidence limits. The violet shading represents 
constraints from the maxBCG cluster catalog, blue those from 
WMAP-5, and the yellow their combination; they have been 
adapted from Allen et al. (2011). The green contours give the to¬ 
mographic peak-count constraints for a Euclid-\i\i& survey cov¬ 
ering 15,000 square degrees. 


errors. These additional studies will quantify the extent to which 
peak counts are complementary to cosmic shear measurements. 

The primary technical challenge in application of peak 
counts is the production of large suites of numerical simulations 
to calculate both the expected mean number of peaks and their 
covariance matrix over the cosmological parameter space. It is 
not, however, unique to the counts: all lensing studies face the 
same challenge because valuable signal, even in the two-point 
statistics of cosmic shear, originates in the non-linear regime. 
We therefore expect peak counts to accompany the more stan¬ 
dard lensing measures in application to large lensing surveys. 
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